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The article describes an experimental method that allows to estimate the inhomogeneous and 
homogeneous linewidths of the photoluminescence band of a point defect in an amorphous solid. 
We performed low temperature time-resolved luminescence measurements on two defects chosen as 
model systems for our analysis: extrinsic Oxygen Deficient Centers (ODC(II)) in amorphous silica 
and ¥^ centers in crystalline Lithium Fluoride. Measurements evidence that only defects embedded 
in the amorphous matrix feature a dependence of the radiative decay lifetime on the emission energy 
and a time dependence of the first moment of the emission band. A theoretical model is developed 
to link these properties to the structural disorder typical of amorphous solids. Specifically, the 
observations on ODC(II) are interpreted by introducing a gaussian statistical distribution of the 
zero phonon line energy position. Comparison with the results obtained on crystalline defects 
strongly confirms the validity of the model. By analyzing experimental data within this frame, we 
obtain separate estimations of the homogenous and inhomogeneous contributions to the measured 
total linewidth of ODC(II), which results to be mostly inhomogeneous. 

PACS numbers: 71.55.Jv, 61.72.jn 78.55.Qr, 78.47.Cd 



I. INTRODUCTION 

The physics of color centers embedded in a solid ma- 
trix is a fundamental and interesting scientific field both 
from the point of view of basic physics and for their wide 
technological applications as modifiers of the macroscopic 
physical properties of solids, (e.g. optical transparency, 
refractive index, electrical resistance, and so on).^'^ Sev- 
eral experimental evidences have led to a general agree- 
ment on the fact that the properties of point defects may 
be significantly different depending on the crystalline or 
amorphous structure of the solid they are embedded iufS. 
Indeed, in a crystal each member of an ensemble of iden- 
tical defects experiences the same local environment. As 
a consequence, every spectroscopical property of the en- 
semble of defects, such as the lineshape of the related 
absorption or photoluminescence (PL) bands, can be in- 
terpreted as a property of the single center, and is re- 
ferred to as homogeneous. The homogeneous absorption 
linewidth is mainly determined by the electron-phonon 
interaction and it is related to other important physi- 
cal properties of the defect, such as the Huangh-Rhys 
factor and the phonon vibrational frequencies.-'- On the 
other hand, defects in an amorphous solid are believed 
to feature site-to-site statistical distributions of the spec- 
troscopic properties due to the disorder of the surround- 
ing matrix. Hence, the lineshapes of their optical bands 
are characterized by an inhomogeneous broadeningj^i^i^ 
which reflects the degree of disorder of the amorphous 
solid and concurs, together with the homogeneous effects, 
to determine the overall spectroscopic signature of the 
color center. 

Many experimental approaches have been proposed to 



estimate the homogeneous and inhomogeneous contribu- 
tions to the experimental linewidth of an optically ac- 
tive center: exciton resonant luminescence, resonant sec- 
ond harmonic scattering, femtosecond photon echo, spec- 
tral hole burning and site-selective spectroscopy.— '^'^'^'^ 
However, the issue is still open since none of these tech- 
niques is applicable to the whole variety of inhomoge- 
neous physical systems of interest. For instance, in amor- 
phous solids site-selective spectroscopy can been success- 
fully applied only to defects which allow the the direct 
observation of the zero phonon line (ZPL) by virtue 
of a weak coupling with the vibrational modes of the 
matrixi^ 



In this paper we propose a new experimental approach 
to this problem, which allows to estimate the inhomo- 
geneous and homogeneous linewidths based on mapping 
the variations of the radiative decay lifetime within an 
inhomogeneously broadened luminescence emission band 
by time-resolved laser-excited luminescence. To this pur- 
pose, in the next section we first describe an adapted ver- 
sion of the theoretical treatment of the optical properties 
of a point defect in a solid, which takes into account the 
effects of heterogeneity in amorphous systems. Next, in 
the experimental section we demonstrate that the pre- 
dictions of our model are consistent with the results of 
measurements performed on two model point defects, one 
in a crystal solid and the other one in a glass. Finally, we 
use the theoretical model to estimate the inhomogeneous 
and homogeneous widths of the two model defects and 
to obtain other physical parameters of interest. 
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II. THEORETICAL DESCRIPTION OF 
OPTICAL DEFECT PROPERTIES 

We briefly review the standard theoretical description 
of the optical properties of a point defect in a crystalr^i^ 
in order to adapt it later to the case of amorphous sys- 
tems. In addition to the crude Born-Oppenheimer and 
Franck- Condon approximations, we suppose the defect 
to be coupled with only one vibrational mode of the 
solid matrix of frequency Wp, assumed to be the same 
for ground and excited electronic states. The frequency 
iOp can be regarded also as the mean frequency of the 
vibrational modes of the solid or can be thought as the 
effective phonon frequency coupled with the electronic 
transitionii^ 

In this frame, the absorption cross section il{E) of a 
defect as a function of the excitation energy E at the 
absolute zero temperature is given hj)^^ 

n{E) =pJ2\Mok\''E ■ 5[E - (Eo + khcjp)] (1) 

k 

where S indicates a Dirac delta function, and the sum- 
mation is carried out over the vibronic transitions link- 
ing the ground electronic state with zero phonons to- 
wards different vibrational sub-levels (fc) of the electronic 
excited state, spaced by hujp. The energy value Eq 
(zero phonon line) is the absorption transition without 
emission or absorption of phonons. Mok is the over- 
lap integral between nuclear wave functions associated 
to the ground and excited states, while (3 is given by: 

/? = i (^^] i^J-lDp where D is the matrix ele- 

ment of the electric dipole operator between the ground 
and excited electronic states, and gi is the degeneracy of 
the lower electronic state. The effective field correction 

accounts for the polarization effect induced 

by the external field on the solid.-^ - We assume here that 
the refraction index n is constant in the electromagnetic 
range investigated. In the harmonic approximation for 
the vibrational sublevels relative to the ground and ex- 
cited states, the |Mofep coefficients are given by a Poisson 
distribution!^ 



-H 



(2) 



where H is the Huangh-Rhys factor, expressing the num- 
ber of phonons emitted by the system after absorption 
of a photon while relaxing to the ground vibrational sub- 
state of the excited electronic levelii Given a population 
of identical defects, the envelop of the 6 functions in Eq. 
(fT|) describes their characteristic homogeneous absorption 
lineshape, with (aside from the effect of the factor E) a 
EAbs = Eo + Hhujp first moment and a — \/Hhujp 
width. 

After relaxation towards the bottom of excited 
electronic state, the system can relax back to the 



ground state by spontaneous photon emission (photo- 
luminescence) . The following relationship of mirror sym- 
metry links the absorption Q.{E) and luminescence L{E) 
band shapesi^ 



L{E) n(2Ea - E) 

(X 



E^ 



1E{) — E 



(3) 



The energy difference between absorption and emission 
peaks {Stokes shift) is linked to the Huangh-Rhys fac- 
tor and results to be 25* = 2Hhu}p. Using the mirror 
symmetry Eq. ([3]) and Eq. ([T]) we obtain: 

L{E) cx \Mok\^E^ ■ S[E - {Eo ~ kUtUp)] (4) 



which represents the homogeneous emission lineshape, 
with (aside from the effect of the factor E^) a E^m = 
Eq — Hhujp first moment and a aho width. Expression 
dH) does not take into account the dependence from the 
excitation energy within the absorption band. This is 
based on experimental results and it will be discussed 
later. 

The PL radiative lifetime r is linked to the absorption 
profile by the Forster's equatiom^^^ 



1/r 



[{2Eo-E)^^dE 
9u J E 



(5) 



where gu is the degeneracy of the upper electronic state. 
Combining Eq. ([5]) and Eq. ^ we obtain the decay rate 
1/r: 



l/T^-lY,\Mok?{Eo-khujpf 



(6) 



where 7 = ^2^2 g^^' "^^^ cubic dependence appearing 
in the above expression is a direct consequence of the 
relation between Einstein coefficients for absorption and 
spontaneous emission, which forms the basis of Forster's 
equation. Eq. ([6]) can be approximated by neglecting the 
contributions far from k ^ H , thus obtaining: 



1/r = -f{Eo ~ Sf 



(7) 



This expression shows that the decay rate is proportional 
to 7 and approximately to the third power of the first 
moment of the emission band. 

Summing up, the global expression for the lumines- 
cence of a population of identical point defects in a solid 
matrix as a function of the spectral position E and time 
t after an exciting light pulse {homogeneous shape) is: 

L{E, t) cx 7^ |A/ofcp£;3e-*/^ • 5[E - {Eo - khcop)] (8) 



This expression assumes that non radiative channels from 
the excited state are absent. As we see from Eq. ([SJ, 
the shape and kinetics of the homogeneous luminescence 
band are completely characterized by four parameters: 



3 



Eq (the ZPL position), hcup (the phonon energy), 7 (pro- 
portional to \D\'^) and H (the Huangh-Rhys factor). H 
and hujp can be expressed in terms of the half Stokes shift 
S and of the homogeneous width aho- ^^p~'^ho/ ^ ^^'^ 
H = S'^/ct^q. In this way, expression ([8]) can be alter- 
natively regarded as depending on the four parameters 
Eq, 5, Uho, 7, thus being indicated by the expression: 
L{E,t\Eo,S, (Thojl)- 

For defects in an amorphous matrix, we can argue the 
hypothesis of a population of identical defects to fail. 
Indeed, each point defect interacts with different envi- 
ronments and it is possible that this conformational het- 
erogeneity causes a site-to-site statistical distribution of 
one or more of the homogeneous properties of single de- 
fects. The simplest model we can put forward to take 
into account the disorder effects is to introduce a gaus- 
sian distribution of the ZPL position i?o, peaked at Eq 
and with an inhomogeneous width (Ti„; in this scheme, 
7, 5*, and auo are still considered as undistributed pa- 
rameters. Within these hypotheses, the global PL signal 
L*{E,t) emitted by the ensemble of non-identical point 
defects can be now expressed as the convolution of the 
homogeneous shape L{E,t) with the inhomogeneous dis- 
tribution of Eq' 



L*{E,t\Eo,ain,S, ahon) « 
L{E, t\Eo, S, aho, 7) • e dEo 



(9) 



Eqs. (Ill) and ([9|) lead us to predict a difference between 
the PL signals of defects in crystalline and amorphous 
solids. Indeed, when the inhomogeneous broadening am 
is almost zero, as expected for point defects in a crys- 
talline matrix, Eq. ([5]) has to be used, and the radiative 
lifetime r should be independent from the spectral posi- 
tion at which it is measured within the emission band. 
In fact, r is expressed by Eq. ([6]), so being a function of 
the homogeneous parameters Eq, 7, S, and aho, which 
are expected to be the same for all defects in the solid. 
In contrast, in an amorphous solid a PL band due to 
an ensemble of point defects can be thought as arising 
from the overlap of several bands with different Eq as 
described by Eq. and thus featuring different life- 
times. Hence, when cri„ is comparable with aho it should 
be possible to experimentally observe a dispersion in t by 
measuring the decay of the PL signal at different emis- 
sion energies. Also, the shape of a band arising from the 
overlap of sub-bands with different lifetimes should vary 
in time, so that the position of its first moment Mi{t), 
calculated by the usual expression: 



Mi{t) 



J E L*{E,t)dE 
J L*{E,t)dE 



(10) 



should depend on time. Therefore, both the dispersion of 
r within the emission band and the time dependence of 
the first moment can be used in principle as experimental 
probes of inhomogeneous effects. 



It is worth noting that according to Eq. ([7]) , r strongly 
depends on the first moment of the emission band, E^.^ = 
Eq — S, and more weakly on 7. This leads to Eq as 
the parameter of choice to be distributed in our model. 
Moreover, a gaussian distribution of Eq was experimen- 
tally demonstrated for the non-bridging oxygen hole cen- 
ter point defect in silica, for which the zero-phonon line 
can be directly observed by site-selective spectroscopy at 
low temperatures^ii^ On the other side, we acknowledge 
that similar predictions can be obtained by introducing 
a distribution of the half Stokes shift S with an undis- 
tributed Eq . Data reported later on in this paper do not 
allow to discriminate between these two possibilities. 

Finally, to get further insight into the meaning of Eq. 
^ it is useful to consider the extreme case in which the 
homogeneous width is so narrow to be negligible with 
respect to the inhomogeneous one. In this case, the 
homogeneous lineshape L{E,t) can be approximated as 
S{E - {Eq ~ S)) ■ e-*/^, with t given by Eq. jT]). By 
substituting in Eq. © we get that: 



L*{E, t) cx e-''-^'* -e' 



(11) 



This expression predicts an exponential decay whose r 
depends cubically from the experimental observation en- 
ergy E within the inhomogeneous band. 

In the intermediate situation of non-negligible homo- 
geneous width, Eq. ([9]) deviates in principle from a sin- 
gle exponential decay, as it contains contributions with 
different values of r. However, we verified that the typ- 
ical values of the parameters which will be used in the 
following to fit experimental data (£^0: ""inj S, aho, 7)7 
correspond to predicted decay curves that always remain 
very close to a single exponential for all practical pur- 
poses. From a theoretical point of view, we can define in 
general t{E) as the time in which L*{E, t) (at a fixed E) 
decreases by a 1/e factor from L*{E, 0). With this defini- 
tion, we can summarize the above considerations as fol- 
lows: the t{E) curve (with E varying within the observed 
emission band) is expected to vary progressively from a 
constant value (for a completely homogeneous system) to 
a cubic dependence (for a completely inhomogeneous sys- 
tem) with increasing inhomogeneous/homogenous ratio. 
To check the validity of our model we have performed 
experimental measurements (described in the following) 
on crystalline and amorphous defects. 



III. MATERIALS AND EXPERIMENTAL 
METHODS 

We chose F-type-centers in lithium fluoride (LiF) and 
Oxygen Deficient Centers of the second type, ODC(II), 
in amorphous silicon dioxide (Si02, or silica) as model 
point defects on which testing our approach. Both cen- 
ters feature broad near-gaussian luminescence bands in 
the ultraviolet (UV) range with close decay lifetime val- 
ues {'^S ns), and they have both been widely studied in 
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literature because of their important technological appli- 
cations. Specifically, LiF is a material traditionally em- 
ployed in the production of high-quality optical elements 
to be used in the infrared, visible, and particularly in the 
ultraviolet spectral regions. F- type-centers in LiF (elec- 
tron trapped in anion vacancies) are the subject of active 
investigation in the areas of color center lasers, radiation 
dosimetry and integrated optics (see Ref. [H and refer- 
ences therein). The study of point defects in silica is a 
fundamental technologic problem as well, because their 
presence compromises the optical and electrical proper- 
ties of glasses in their wide uses as optical components, as 
insulators in MOS transistors, and for guiding or process- 
ing light signals (optical fibers and Bragg gratings) 
ODC(II) is a peculiar defect of the amorphous phase of 
Si02,'^'^* thus being an interesting model system to in- 
vestigate the characteristic properties of defects in dis- 
ordered materials with respect to crystalline ones. Its 
microscopic structure consists in an atom bonded to two 
oxygen atoms of the matrix (=X" ) , where X is an atom 
belonging to the isoelectronic series Si, Ge, Sn^iiiii^ Pre- 
vious studies have suggested that the spectroscopic prop- 
erties of ODC(II) are significantly conditioned by inho- 
mogeneous effects J^iiiii^'i^ 

We report measurements performed on two samples: 
the first one is a crystalline Lithium Fluoride sample, 
hereafter denoted as LiF. Prior to any measurement this 
specimen, 5 x 5 x 1.25 mm"^ sized, was irradiated at room 
temperature with electrons of 3 MeV energy, for a to- 
tal dose of 1.5 • 10^ rad. The purpose of irradiation 
was to induce in the sample the formation of lumines- 
cent F-type centers. The second sample is a fused sil- 
ica (commercial name: InfrasilSOl, provided by Heraeus 
Quartzglasr^ and 5x5x1 mm'^ sized), hereafter named 
1301, manufactured by fusion and quenching of natural 
quartz, with typical concentration of impurities of ^^20 
ppm in weighti^S In particular, as-grown 1301 contains 
a ~1 ppm concentration of Ge impurities, due to con- 
tamination of the quartz from which the material was 
produced. Previous studies demonstrated that in the as- 
grown material most of the Ge impurities are arranged as 
Ge-ODC(II) defects (=Ge"); moreover comparison with 
sol-gel silica samples doped with Ge atoms ensures us 
that the contribution to PL of intrinsic ODC(II) defects 
in 1301 sample is negligible J^i^i The optical activity of 
Ge-ODC(II) at low temperature (<100 K) consists in an 
absorption band centered at ~ 5.1 eV which excites a 
fast (lifetime in the ns range) emission band centered at 
~ 4.3 eV, due to the inverse transitiouFiSiiiiS^ 

PL measurements were done in a standard back- 
scattering geometry, under excitation by a pulsed laser 
(Vibrant OPOTEK: pulsewidth of 5 ns, repetition rate 
of 10 Hz, energy density per pulse of 0.30±0.02 mJ/cm^) 
tunable in the UV- Visible range. The luminescence emit- 
ted by the sample was dispersed by a spectrograph (Spec- 
traPro 23001, PI Acton, 300 mm focal length) equipped 
with three different gratings, and detected by an air- 
cooled intensified CCD (Charge-Coupled Device PIMAX, 



PI Acton). The detection system can be triggered in or- 
der to acquire the emitted light only in a given temporal 
window defined by its width [iw] and by its delay to 
from the end of the laser pulse. All measurements re- 
ported here were performed on samples kept at 25 K in 
high vacuum (~ 10~^ mbar) within a He flow cryostat 
(Optistat CF-V, OXFORD Inst.). AU luminescence sig- 
nals in 1301 were acquired with a 300 grooves/mm grating 
with a 2 nm bandwidth, while the signals in LiF were 
measured with a a 150 grooves/mm grating with a 2.5 
nm bandwidth. All the spectra were corrected for the 
spectral response and for the dispersion of the detection 
system. 



IV. EXPERIMENTAL RESULTS 

In Fig. [T]-(a) we show a typical time-resolved mea- 
surement of the PL activity of Ge-ODC(II) in the 1301 
sample, performed at 25 K under laser excitation at 240 
nm (5.17 eV). The PL decay was analyzed by performing 
60 acquisitions with the same integration time tvK=l ns 
but at different delays t o , going from to 60 ns from the 
laser pulse. Fig. [T]-(b) shows the normalized spectra of 
panel (a) in a contour plot and evidences that the first 
moment of the band (continuous line) varies in time. 

In Fig. [l]-(a) we report the signal acquired for tD=0, 
corresponding to the first spectrum in Fig. [l]-(a). The 
PL band of Gc-ODC(II), as acquired immediately after 
the end of the laser pulse, is peaked at ^4.4 cV and has a 
-0.45 eV width (FuU width at Half Maximum, FWHM) 
consistent with literature dataj^^ Completely analogous 
time-resolved measurements were carried out on the PL 
activity of F-type centers in the LiF sample. 

This specimen was excited at 450 nm (2.76 eV) and its 
luminescence was collected by varying Id from to 100 
ns with tiy = l ns. We report in Fig. [ll-(b) the lumines- 
cence signal detected in LiF at io—^- It is apparent that 
the PL signal of LiF comprises two contributions peaked 
at —2.3 eV and —1.8 eV. These signals are known to be 
associated to two different defects, the and F2 cen- 
ters respectively, both consisting in aggregates of F-type 
centersii^i^ In particular, the main —2.3 eV band with 
a —0.27 eV FWHM is due to , consisting in two elec- 
trons localized on three adjacent anion vacanciesJ^ For 
each activity (Ge-ODC(II) and F3+), one can extract the 
time dependence of the first moment of the luminescent 
bands from the time-resolved measurements (e.g. those 
in Fig. [l]in the case of Ge-ODC(II)). Data so-obtained 
are reported in Fig. [5]-(a). The origin of the time scale 
corresponds to t^i^O. We observe that the PL activity in 
silica shows an approximately linear decrease of the first 
moment as a function of time, while this decrease is not 
observed in LiF, where the first moment of the F^ cen- 
ters band has a constant value within experimental sen- 
sitivity. As already discussed in the theoretical section, 
the progressive shift of the PL peak position observed 
for ODC defects can be alternatively understood as a de- 
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FIG. 1; (a) Decay of the luminescence band of Ge-ODC(II) centers in the 1301 sample, excited at 240 nm at T=25 K. (b) 
Normalized data of panel (a) in a contour plot. The continuous line corresponds to the position of the first moment of the PL 
band as a function of time. The position of the first moment at tD=0 is reported (dashed line) as a reference. 
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FIG. 2; Low temperature (25 K) luminescence of Ge-ODC(II) 
in the 1301 sample (a), and of F-centers in the LiF sample (b). 
Both PL bands are obtained by exciting at the maximum 
of the respective absorption bands and acquired for to = 
and with tw~l ns. The continuous line is the result of the 
fitting procedure by our theoretical model. The Poissonian 
homogeneous shape is also shown (see discussion). 



pendence of the luminescence lifetime from the spectral 
position within the emission band. Hence, in Fig. [3]-(b) 
we report the values t{E) of the PL hfetime as a func- 
tion of the emission energy. The lifetimes were estimated 
for both PL activities by least-square fitting data from 
time-resolved spectra ([I]- (a)) at different emission ener- 
gies with an exponential function {I{t) = /(O)e^*/'^).'^'^ 
At this temperature (25 K), the decays are purely expo- 
nential for both activities.— Fig. [3]-(b) shows that the 
lifetime of Ge-ODC(II) centers in silica strongly varies 
within the emission band: r goes from ^--^7.0 ns to ~10.7 



ns. A similar behavior for Ge-ODC(II) was observed also 
under excitation by synchrotron radiation. On the con- 
trary, the lifetime of F^^ centers is almost constant in the 
observed range of emission energies. The above results 
were obtained exciting at the absorption peak for both 
PL activities. Although we performed the same measure- 
ments for different excitation energies within the absorp- 
tion band, only a very weak dependence from this param- 
eter was evidenced, consistently with previous results 



V. DISCUSSION 

The results in Fig. [3] qualitatively confirm the predic- 
tions of our theoretical analysis, i.e. that the dependence 
of the lifetime on the emission energy or, equivalently, the 
progressive red-shift of the emission peak with time, are 
characteristic features of luminescent defects embedded 
in a glassy matrix, as opposed to "crystalline" defects. 
We stress that the non-radiative decay channels are al- 
most completely quenched for both PL signals at the 
temperature at which the experiments were performed 
(25 K) J^i^ As a consequence, it is a very good approx- 
imation to consider the luminescence decay to be purely 
radiative. The main point of the following discussion is 
to fit all experimental data by our model and extract the 
values of the homogeneous and inhomogeneous widths 
of the PL emission bands and other interesting physical 
parameters. 

For both investigated PL activities we have performed 
numerical integration of Eq. ^ to obtain a set of three 
theoretical curves which simultaneously fit i) the shape of 
the PL band at tD=0, ii) the time dependence of the first 
moment (calculated by Eq. pO)) ) and iii) the dependence 
of r on emission energy?^. To increase the reliability of 
the fit procedure, the half Stokes shift S was fixed to 
the value obtained experimentally by measuring the dif- 
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FIG. 3: (a) First moment of the emission band of Ge-ODC(II) (white circles) and of (white squares), (b) Decay hfetime 
as estimated by fitting with an exponential function data at different emission energies within the emission band of the Ge- 
ODC(II) centers in the 1301 sample (white circles) and of the Fj^ in the LiF sample (white squares). The continuous lines are 
the results of the fitting procedure by our theoretical model (see discussion). 



ference between the spectral positions of the absorption 
and emission peaks: S'=0.38 eV and 3=0.24: eV in silica 
and LiF respectively. In this way, the fitting procedure 
was performed by varying only four free parameters, Eq, 
o'in, cTho, 7- From the experimental point of view, the 
vibrational sub-structure of homogeneous luminescence 
bands cannot usually be resolved due to the bandwidth 
of the measuring system and to further broadening ef- 
fects due for instance to the coupling with several low 
energy modes. To take into account this effect, the ho- 
mogeneous lineshape, Eq. ([5]), was convoluted with a 
gaussian distribution of a narrow width Tiujp before being 
inserted into Eq. ([5]). 

The continuous lines in Fig. [5] and [3] represent the 
results of our fitting procedure. It is worth underlin- 
ing the goodness of the fit, obtained using only four pa- 
rameters, and considering especially that data in Fig. 
[3] take into account simultaneously all data acquired in 
a time-resolved PL measurement (typically ^600 spec- 
tral positions for each of the ^^100 temporal acquisi- 
tions of Fig. [T|) . Table U summarizes the best param- 
eters obtained via our fitting procedure for the two in- 
vestigated PL activities. From data in Table [T] we can 
also calculate the Huang-Rhys factor H = S'^/ct^q, the 
vibrational frequency hujp = af^^jS, the total width 
(from a^g^ = tT^„ + cr^^jo),"^^ and finally the parameter 
A = cr infer tot which estimates the degree of inhomogene- 
ity. All these quantities are reported in Table Has well. 
As expected, A is very small for the LiF defects in com- 
parison with the amorphous ones: ~ 3% against ^ 78%. 
These values correspond to am being about 0.2 times and 
2 times aho, in LiF and Si02 respectively. We note that 
the inhomogcncous broadening in the crystalline sam- 
ple is not exactly zero; beside the approximations in our 
model, we note that a real crystal is always distorted 
by some dislocations, strains or other imperfections dis- 
tributed at random into the matrix. The obtained value 
of A for Ge-ODC(II) shows that for a defect embedded in 



a glassy matrix the inhomogeneous width can be promi- 
nent with respect to the homogeneous one. This conclu- 
sion may be at variance with previous suggestions that 
aho and am are typically comparableJi^ 

In Fig. [5] we also show the discrete Poissonian homo- 
geneous lineshape of width a ho and ZPL position Eq, as 
obtained by our fit procedure for both investigated ac- 
tivities. As already pointed out, the crystalline PL band 
is completely described by the homogeneous shape,— 
whereas the silica PL band is not reproduced without 
taking into account inhomogeneous effects. It is also 
worth noting that the value hu!p=23±6 meV obtained via 
our fitting procedure is very close to the value of 26 ± 2 
meV found for the same defect by the analysis of the 
temperature dependence of the experimental absorption 
linewidthi^ Moreover, hujp is in good agreement with 
experimental and computational works on silica glasses 
which predict the presence of vibrational modes of low 
frequencyj ^^i^^ These agreements further confirm the cor- 
rectness of the present value of aho and consequently of 
our analysis. 

To show the accuracy of our fitting procedure in deter- 
mining A, in Fig. [3]we compare the experimental lifetimes 
of Ge-ODC(II) in the 1301 sample with the predictions 
of our model obtained for different A values. The theo- 
retical t{E) curves are obtained by keeping atot fixed to 
the value which best fits the overall experimental shape 
of the PL band. This analysis clearly evidences a contin- 
uous transition from constant lifetimes for A=0 (that is 
a completely homogeneous PL band) , to an inverse cubic 
dependence of t from emission energy for A=l (that is 
a completely inhomogeneous PL band) , as anticipated in 
the theoretical section. Finally, the oscillator strength / 
reported in Table U is calculated using)^ 



where nie and e are respectively the mass and the charge 
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TABLE I: Upper section: best fitting parameters obtained by our theoretical model for the investigated PL activities. Lower 
section: Values of A, atot, hijjp, H , and /, as calculated from best fitting parameters 





[eV] 


cTin [meV] 


aho [meV] 


S[eV] 


7 [10*' eV-^s-''] 


1301 


4.70±0.05 


177±10 


93±12 


0.38±0.02 


1.41±0.09 


LiF 


2.50±0.02 


20±10 


109±6 


0.24±0.02 


10.0±0.6 




A(%) 


atot [meV] 


hujp [meV] 


H 


/ 


1301 


78±5 


200±10 


23±6 


17±5 


0.073±0.010 


LiF 


3±2 


111±6 


51±7 


5±1 


0.32±0.04 




3.8 4.0 4.2 4.4 4.6 4.8 
Energy (eV) 



FIG. 4: Experimental decay lifetime at different emission en- 
ergies for Ge-ODC(II) point defects (white circles). Lifetime 
as predicted by our model for different values of the parameter 
A (continuous lines) . Dashed line represents the extreme case 
of 1/-B^ dependence (see discussion) . The arrow indicates the 
direction of increasing A. 

of electron. We have substituted in Eq. the value 

of IDp, calculated from the fitting parameter 7, and we 
have used for EAbs the value Eq + S. In regard to the 

effective field correction, the term — ( ] calculated 

within the Onsagcr modelf^i^ results to be close to unity 
both in Si02 (n~1.5) and in LiF (n~1.4) in the investi- 
gated spectral range. The oscillator strength found here 
for Ge-ODC(II) in silica is consistent with the range of 
values reported in UteratureiJ^ 0.03-0.07. For centers 
in LiF our result is close to 0.2 reported in ref.l29l. 

The main assumption of our model that all amorphous 
effects can be completely accounted for by a simply gaus- 
sian distribution of a single homogenous parameter (i.e 
zero phonon line) is strongly corroborated by the excel- 
lent agreement between theoretical curves and data. On 
the other side, a distribution of the emission peak Eq — S 
is strongly suggested a priori by the almost Einstein- 
like proportionality of 1 /r on E"^ shown by experimental 
data in Fig. [H Finally, it is important to note that in this 
scheme 7 and thus |Z3p are assumed as undistributed pa- 



rameters. This means that the oscillator strength given 
by Eq. can be distributed only as a consequence 

of the variations of Eai,s associated to different homoge- 
neous absorption sub-bands. 



VI. CONCLUSIONS 

We have investigated the inhomogeneous properties of 
point defects in a glassy matrix via mapping by time- 
resolved PL the dependence of the radiative decay life- 
time on emission energy. We propose a theoretical model, 
based on an extension of the standard theory of the opti- 
cal properties of point defects, incorporating a statistical 
distribution of the zero phonon line to account for the ef- 
fects of the non-equivalent environments probed by each 
point defects in an amorphous matrix as opposed to a 
crystalline one. This model enlightens a direct connec- 
tion between the dispersion of the radiative decay lifetime 
within a luminescence band as a function of emission en- 
ergy and the inhomogeneous properties of defects in a 
glassy environment. To confirm our prediction we have 
experimentally studied the luminescence of Oxygen De- 
ficient Centers in silica and of aggregates of F-centers 
in a crystalline sample of LiF. The model is able to fit 
all experimental data and to provide an estimate of the 
ratio A = crf^/a^^^ between the inhomogeneous and the 
total width, namely -78% for ODCs and -3% for F+. 
Finally, our model allowed us to determine the homoge- 
neous parameters of ODC and Fj^ centers: homogeneous 
width, oscillator strength, Huangh-Rhys factor and the 
frequency of the vibrational local mode. 
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